48?^) 


IttViSRnON  EQU1*«NT  NlMUIRRCTUREO  B¥  THE  MAEKR^ 
COIVORATIOR,  1^:- 


P A C KemMMty  and  T A Jones 


ADMIRALTY  SURFACE 
PORTSMOUTH  HARTS 


ESTABUSWieKnF 


’¥'ij 


Copyright  (DControlleMilSO,  London  197^ 


‘DTSTRtwnTIQN  STATEMENT  t 
Apptovwl  lof  public  wUoMt 


To  all  Recipients 


ASWE 

Portsdown,  Cosham 
Portsmouth,  Hants 
P06  4AA 

XMS5/INF04.1 

September  1976 


ASWL  TECHNICAL  REPORT  TR  76005  (UNLIMITED) 

THE  STATIC  PERFORMANCE  OF  THE  MX902  SATELLITE 
NAVIGATION  EQUIPMENT  MANUFACTURED  BV  THE  MAGNAVOX 
CORPORATION,  USA 

ERRATA 


PAGE  (iii)  APPENDIX  K Page  Nos  K1 , K2  should  read  K1-K4. 

PAGE  1 , PARAGRAPH  1 Second  line  should  read 

" available  to  the  Admiralty  Compass  Observatory". 

PAGE  4,  LINES  26  AND  28  " are  based  on  the  rms " should  read 

" is  based  on  the  rms 

PAGE  6,  PARAGRAPH  5 Sixth  line  should  read  "deposited  it  in  a file  SATR". 

PAGE  7,  PARAGRAPH  5.1.1  Last  paragraph.  Line  2 

"distributon"  should  read  "distribution". 

PAGE  9,  LINE  3 "x  is  the  latitude"  should  read 

"x.  is  the  latitude". 

J 

PAGE  9,  LINE  4 "y  is  the  latitude"  should  read 
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ABSTRACT 

^he  report  contains  details  of  an  investigation  into  the  accuracy  and  error 
patterns  obtained  with  the  satellite  navigation  equipment,  MX902,  produced  by 
the  Magnavox  company  of  the  United  States.  The  readings  were  statistically 
analysed  and  found  to  be  in  good  agreement  with  the  manufacturer's  predicted 
performance.  The  readings  were  also  studied  with  a view  to  finding  a 
relationship  between  positional  error  and  any  of  the  other  calculated  variables 
such  as  offset  frequency  and  number  of  iterations;  no  such  relationship  was 
found. 
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1.  INTRODUCTION 


The  satellite  navigation  equipment,  MX902,  manufactured  by  the  Magnavox 
Company,  was  made  available  at  the  Admiralty  Compass  Observatory  for  a short  time 
before  being  fitted  on  board  a Royal  Naval  vessel.  This  opportunity  was  taken  to 
assess  the  static  accuracy  and  error  patterns  associated  with  the  equipment.  No 
mobile  trials  were  conducted,  thus  the  effect  of  vehicle  motion  could  not  be 
determined.  The  fixes  obtained  were  analysed  by  means  of  a set  of  computer  programs 
written  in  standard  Fortran  IV,  and  these  were  run  on  Comshare's  Rank  Xerox  RXDS 
Sigma  9 computer.  Programs  from  Comshare's  statistical  package  (TACTICS)  were  also 
used  in  the  analysis. 

The  assessment  presented  here  is  an  extension  of  the  work  documented  in  ASWE 
Technical  Report  TR-74-19  (Ref  1).  Since  that  report  was  written  more  readily 
accessible  computing  facilities  have  been  made  available,  via  a contract  with  the 
time  sharing  bureau,  and  this  has  necessitated  the  re-writing  of  the  analysis 
programs  used  initially.  Communication  was  via  two  GPO  modem  links,  one  to  a 
10  character  per  second  teletype  and  the  second  to  a 30  character  per  second 
Tektronic  4010-1  visual  display  unit. 

Comshare  have  a batch  facility  on  their  computer,  and  using  this  to  run  the  larger 
analysis  loads  overnight  gave  a considerable  saving  in  cost. 

Three  types  of  position  fix  were  available  for  analysis:- 

(a)  Position  fixes  which  the  system,  basing  its  decision  on 
predetermined  parameters,  considered  adequate  to  use  as  a valid 
"position  up-date". 

(b)  Position  fixes  which  the  system,  using  the  same  parameter  as  in 
(a),  considered  poor  and  not  to  be  relied  upon. 

(c)  "Velocity  North"  fixes  which  were  calculated  as  a result  of  a 
re-estimate  of  the  assumed  northerly  velocity  of  the  receiving 
aerial,  based  on  the  fix  provided  in  (a)  above.  Since  the  receiver 
used  in  the  assessment  was  sited  in  a laboratory,  V|^  was  initially 

assumed  to  be  zero. 

In  practice  other  factors,  apart  from  an  actual  error  in  assessing  or  measuring  V^^, 
can  give  the  same  result  as  an  error  in  V,,.  A correction  to  this  value  can 
therefore,  in  theory,  improve  the  quality  of  the  resulting  position  fix. 

2.  OUTLINE  OF  THE  NAVY  NAVIGATION  SATELLITE  SYSTEM 


The  United  States  Navy  Satellite  System,  known  as  TRANSIT,  was  developed  by 
the  US  Navy  in  conjunction  with  Johns  Hopkins  University  during  the  early  1960s, 
and  became  operational  in  1964.  At  the  time  of  this  assessment,  the  system 
consisted  of  six  operational  satellites  in  near  polar  orbit  and  with  an  average 
altitude  of  1000  km.  Each  satellite  therefore  circumnavigates  the  earth  once  every 
108  minutes.  For  an  observer  positioned  at  either  pole,  a satellite  would  appear 
over  the  horizon  on  average  every  18  minutes,  as  each  orbit  would  pass  almost 
directly  overhead.  As  the  observer  moves  towards  the  equator,  the  frequency  at 
which  the  satellites  appear  over  the  horizon  decreases  as  a substantial  number  of 
orbits  are  masked  by  the  earth,  and  do  not  rise  over  the  local  horizon.  At  the 
equator,  the  average  period  between  observable  satellite  passes  is  approximately 
90  minutes.  As  the  satellites  are  launched  they  are  given  an  orbit  which  is 
displaced  inlongitude  to  the  others,  and  ideally  these  angular  separations  should 
be  equal.  As  more  satellites  are  introduced  to  fill  in  the  gaps  the  symmetry  must 
be  upset,  but  this  would  occur  with  time  in  any  case.  The  altitude  of  each 
satellite  differs  slightly  and  this  will  cause  them  to  have  differing  precession 
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10.  CONCLUSIONS 


rates  and  orbital  periods.  This  means  that  with  time,  the  orbital  separations  will 
slowly  change  and  give  rise  to  a slow  beat  producing  what  is  known  as  bunching. 

The  effect  of  bunching  is  particularly  pronounced  ^t  equatorial  latitudes,  where 
for  certain  locations  it  can  cause  periods  of  up  to  seven  hours  when  no  satellite 
passes  are  observable.  This  period  can  then  be  followed  by  a succession  of 
observable  passes  over  a fairly  short  interval.  After  a few  days  this  effect  will 
disappear,  as  the  orbits  begin  to  spread  out  again. 

Each  satellite  is  provided  with  a "Digital  Store",  a stable  oscillator,  an  S-Band 
receiver,  two  low  power  transmitters  and  an  antenna  system.  The  satellites  are 
also  gravity  stabilised  to  ensure  that  the  antenna  systems  are  always  directed 
towards  the  earth.  As  the  satellites  orbit  the  earth,  they  continuously  transmit 
two  stable  frequencies,  which  are  phase  modulated  with  data  stored  in  the 
satellite's  memory.  This  data  describes  the  actual  satellite  position  at  the  start 
of  each  "Universal  Time"  two  minute  period.  One  complete  set  of  data  takes  just 
less  than  two  minutes  to  transmit,  and  is  known  as  the  "Satellite  Message".  The 
transmissions  from  each  satellite  are  monitored  in  various  parts  of  the  world  by 
ground  stations,  and  these  are  equipped  to  determine  small  variations  between  the 
actual  orbit  and  that  described  by  the  Satellite  Message.  Because  of  the  satellite's 
environment  small  orbital  changes  do  take  place,  and  a master  ground  station 
equipped  with  a suitable  transmitter  is  able  to  refresh  the  variable  data  in  the 
satellite's  digital  store  by  a remote  injection  procedure.  In  practice  this 
updating  occurs  every  12  to  14  hours  for  each  satellite. 

The  frequency  standards  used  in  the  satellites  are  good  quality  crystal  units,  and 
are  suitably  synthesised  to  provide  both  the  digital  timing  waveforms  and  the 
transmitter  frequencies.  The  frequencies  transmitted  by  the  satellite  are  nominally 
400  MHz  minus  32  kHz,  and  150  MHz  minus  12  kHz.  The  difference  between  400  MHz  and 
the  transmitted  frequency  (minus  32  kHz)  is  known  as  the  OFF-SET  FREQUENCY  and  is  a 
quantity  determined  by  the  receiving  equipment.  In  most  receiving  sets  it  is 
displayed  in  cycles/minute  (32  kHz  = 1920000  c/min)  and  is  initially  assumed  to  be 
unknown.  Although  the  higher  transmitted  frequency  is  nominally  399.968  MHz,  every 
satellite  will  vary  slightly  as  the  frequency  standard  experiences  thermal  cycling 
and  aging.  The  local  standard  in  the  receiver,  against  which  the  satellite 
frequency  is  measured,  will  also  drift  and  therefore  must  be  accounted  for  in  the 
computation.  The  aging  of  the  satellite's  frequency  standard  will  normally  be  slow, 
and  therefore  given  sufficient  passes  it  should  be  possible  to  identify  frequency 
bands  which  move  slowly  with  time.  In  the  MX902-3  systems,  the  frequency  offset  is 
calculated  and  used  to  determine  the  deviation  of  the  receiver  local  oscillator  from 
its  design  value.  This  is  then  provided  as  N parts  in  10^,  and  was  evaluated  in  the 
following  analysis  as  a possible  position  error  monitor. 

If  the  navigator  is  to  use  the  information  available  from  the  satellite  to  obtain  a 
position  fix,  a scheme  has  to  be  devised  to  relate  the  satellite's  position  (which 
is  known)  to  that  of  the  observer.  This  was  achieved  by  developing  special 
receiving  equipment  that  would  measure  the  change  in  range  of  an  observable 
satellite,  with  respect  to  a mobile  ground-based  unit,  as  the  satellite  orbits  the 
earth.  If  the  observer  remains  stationary,  the  change  in  range  would  be  due  only  to 
the  satellite  motion.  If  the  observer  moves,  a component  of  the  range  change  would 
be  due  to  his  motion  over  the  ground  and  must  be  taken  into  account.  The  system 
used  in  the  early  receivers  was  to  measure  the  change  in  range  over  each  two  minute 
satellite  period,  which  is  well  defined  by  timing  marks  at  the  beginning  of  each 
period.  This  change  in  range  over  the  two  minute  period  was  measured  by  comparing 
the  received  satellite  frequency  with  a local  reference,  and  determining  the  number 
of  "doppler"  cycles  occurring  in  the  measurement  period.  By  this  means  a direct 
measurement  of  the  change  in  range  can  be  made.  To  obtain  a fix  three  range  changes 
are  required  in  order  to  determine  three  unknowns  - namely  LATITUDE,  LONGITUDE  and 
OFF-SET  FREQUENCY. 
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the  VN  solution  could  then  be  determined  in  a dynamic  environment.  It  is  possible 


The  disadvantage  with  the  system  as  described,  is  that  any  loss  of  signal  during  the 
two  minute  period  would  upset  the  measurement  for  the  whole  period.  To  overcome 
this,  the  current  equipments  further  divide  the  two  minute  periods  into  a number  of 
shorter  periods.  In  the  MX902  equipment,  the  two  minute  period  is  divided  into  five 
periods:  four  of  23  seconds  and  one  of  27.6  seconds.  These  in  turn  are  built  up 
from  the  basic  counting  interval  used  by  the  receiver  of  ^.6  seconds.  Since  these 
basic  counting  periods  are  available  to  the  receiver  there  is  great  flexibility  in 
the  way  they  are  compiled,  and  special  schemes  could  easily  be  devised  to  limit  the 
effect  of  wavewash  in  submarine  installations  or  for  an  obstructed  aerial. 

The  above  scheme  works  because,  within  the  area  of  interest,  the  measured  range 
changes  are  peculiar  to  only  one  location.  By  successive  iterations  from  an  assumed 
position  the  actual  position  may  be  resolved.  One  measure  of  the  quality  of  the  fix 
are  the  residuals,  which  are  the  rms  values  of  the  differences  between  the  measured 
and  finally  calculated  range  differences  used  in  the  fix  computation. 

One  novel  feature  of  the  MX902  series  of  equipments,  is  a program  technique  known 
as  the  "Velocity  North"  solution.  The  requirement  to  know  the  vehicle's  speed  over 
the  ground  has  already  been  briefly  discussed,  and  any  error  in  the  assessment  of 
this  quantity  will  give  rise  to  erroneous  position  fixes.  The  position  error 
patterns  are  quite  systematic  with  velocity  errors,  but  the  problem  is  that  the 
actual  velocity  errors  are  not  known.  If  they  were,  correction  to  the  velocity 
inputs  could  be  made  to  prevent  the  position  errors  arising  in  the  first  place. 
Typically  an  error  of  1 knot  in  the  Velocity  North  component  of  vehicle  motion  will 
give  rise  to  a position  error  in  longitude  of  0.2  nm,  and  an  error  of  1 knot  in  the 
Westerly  Velocity  will  give  an  error  of  0.05  nm  in  latitude.  To  overcome  this 
problem,  Magnavox  have  devised  a technique  of  using  the  normal  position  fix  to 
assess  the  actual  vehicle  speed,  and  then  use  this  to  recompute  a further  position 
fix.  All  the  parameters  that  could  give  rise  to  Velocity  North  type  errors,  are 
assumed  to  be  due  to  an  error  in  the  V^^  estimate.  It  was  postulated  that  this  may 

cause  problems,  and  the  assessment  presented  in  this  report  addresses  this.  The 
"trial  receiver"  was  stationary  throughout  the  data  gathering  periods,  and  therefore 
no  V|^  errors  should  have  occurred.  Any  non-velocity-generated  errors  treated  as 

velocity  errors  would  therefore  be  assumed  to  be  systematic,  and  they  may  well  in 
fact  be  random. 

3.  DATA  PREPARATION 


The  following  section  details  the  methods  used  to  convert  the  available  data 
into  a form  suitable  for  use  by  the  statistical  programs. 

3.1  Initial  Preparation 

The  initial  data  is  in  the  form  of  a teletype  log,  see  Figure  1.  Each  fix 
record  consists  of  the  title  of  the  program:  "400  Fix  ACO-SY-75030"  in  the  case  of 
a standard  fix,  and  "400  VN  Fix  ACO-SY-75030"  in  the  case  of  a Velocity  North  fix. 

The  second  line  is  a record  of  the  Doppler  periods,  each  two  minute  Doppler  period 
being  split  into  five  sub-periods,  four  of  23  seconds  and  one  of  28  seconds.  Each 
character  in  a line  represents  one  sub-period,  and  they  have  the  following  meanings:- 

1 Good  data  received. 

0 Very  bad  or  no  data  received. 

? Data  taken  when  satellite  below  7J°. 

= Data  had  some  errors  but  the  computer  was  able  to  edit  them  out. 
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The  line  below  this  contains  two  arrows  (tt)  which  indicate  where  the  centre  of  the 
pass  occurred. 


The  next  two  lines  list  the  various  measured  and  calculated  quantities,  with  their 
values  on  the  line  beneath. 

The  quantities  are:- 

First  line  - DAY  - Julian  day  numoer. 

TIME  - Hours,  minutes,  seconds  in  GMT. 

LAT  - Fix  Latitude  in  degrees,  minutes  and  decimals  of  a minute 
followed  by  the  direction,  N or  S. 

LON  - Fix  Longitude  in  degrees,  minutes  and  decimals  of  a minute, 
followed  by  the  direction,  E or  W. 

ANT  - The  height  of  the  aerial  above  the  reference  geoid  in 
metres . 

HOG  - Ship's  heading  in  degrees. 

SPD  - Ship's  speed  in  knots  (usually  speed  through  the  water). 

NVEL  - Calculated  Velocity  in  a Northerly  direction  in  knots. 

Second  line  - ITER  - The  number  of  iterations  taken  by  the  program  to  calculate 

the  resulting  fix. 

ELEV  - The  maximum  elevation  of  the  satellite  above  the  horizon. 

GEOM  - The  direction  of  the  satellite,  eg  N-W  indicates  the 

satellite  is  travelling  North,  and  is  to  the  West  of  the 
observer. 

SAT  - The  satellite  identifier.  There  are  six  satellites,  30120, 
30130,  30140,  30180,  30190,  30200. 

S-LA  - The  standard  deviation  of  the  latitude  in  minutes  of  arc 
are  based  on  the  rms  doppler  residuals. 

S-LO  - The  standard  deviation  of  the  longitude  in  minutes  of  arc 
are  based  on  the  rms  doppler  residuals. 

S-VN  - The  standard  deviation  of  the  Northerly  Velocity. 

RMS  - Calculated  rms  error  of  fix,  in  nautical  miles,  using  S-LA 
and  S-LO. 

FREQ  - The  offset  frequency  of  the  local  oscillator  with  respect 
to  its  design  value  in  parts  in  10^.  This  should  be  small 
and  not  change  by  more  than  3 units  per  day.  If  a satellite 
pass  results  in  a value  which  differs  by  more  than  a few 
counts  from  the  results  of  other  passes,  then  it  was 
considered  to  be  a poor  position  fix. 

The  S-UPDATE  line  gives  the  differences  in  nautical  miles  between  the  present  fix 
and  the  dead  reckoned  position  fix,  in  latitude,  longitude  and  radial  distance. 

The  C-WHDG  and  C-WSPD  represent  the  calculated  water  heading  and  calculated  water 
speed.  The  calculations  are  based  on  the  drift,  calculated  from  the  S-UPDATE  data, 
divided  by  the  time  since  last  update. 

The  last  line  indicates  whether  the  system  considered  the  fix  sufficiently  accurate 
to  be  used  as  an  update,  this  shown  by  the  word  UPDATE,  or  N-UPDATE  in  the  case  of  a 
poor  fix. 
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The  data  was  punched  onto  paper  tape  which  was  then  fed  into  the  bureau  computer's 
store,  here  it  was  stored  in  a file  called  SATOAT. 

3.2  Wild  Point  Editing 


Since  the  data  was  produced  while  the  receiver  was  at  a fixed  position,  ideally 
the  latitude  and  longitude  should  have  been  invariant,  but  to  prevent  the  analysis 
of  the  data  being  corrupted  oy  "gross  errors",  it  was  considered  necessary  to  edit. 

There  were  two  stages  to  the  editing  procedure,  the  first  perfornied  by  the  analysing 
computer  and  the  second  by  the  investigator.  This  second  editing  takes  place  at  a 
later  stage,  see  Section  4.1. 

3.2.1  Computer  Editing 

The  program  SAT  1 read  the  records  from  SATDAT,  and  transferred  part  of 
this  to  the  file  SATD.  It  also  checked  that  the  number  of  degrees  of 
latitude  and  longitude  were  correct  (this  was  justified  as  the  latitude 
and  longitude  were  approximately  in  the  centre  of  two  integer  degree 
poi nts ) . 

In  the  editing  procedure  two  fixes  were  edited  out,  these  being  the 
standard  and  Velocity  North  fix  from  the  same  satellite  pass. 

Examination  of  this  fix  showed  that  the  local  oscillator  offset 
frequency  in  both  cases  was  very  large,  -3187.6  and  -3044.5  respectively, 
which  indicates  that  the  receivers  had  not  fully  warmed  up.  This  had 
caused  very  large  errors  in  latitude  and  longitude. 

Appendix  A contains  a listing  of  SAT  1,  the  print  out  from  a run  of 
SAT  1 and  a print  out  of  the  two  erroneous  records.  Appendix  B contains 
a listing  of  file  SATD. 

4.  PRIMARY  CALCULATIONS 

Using  the  TACTICS  statistical  programs,  the  data  on  file  SATD  was  first  split 
into  three  partitions.  The  first  contains  those  ordinary  fixes  that  the  computer 
used  to  update  the  current  position,  and  was  called  partition  FIX. UPDATE.  The 
second  contains  ordinary  fixes  that  the  computer  did  not  use  to  update  position,  and 
was  called  partition  FIX. NOUPDATE.  The  third  contains  all  Velocity  North  fixes  and 
was  called  partition  VNFIX.  All  partitions  were  treated  independently. 

The  true  latitude  and  longitude  of  the  receiver  site  were  51°  29.45'  North  and 
0°  33.583'  West,  these  were  determined  by  an  accurate  survey  and  are  based  on  the 
UK  Ordnance  Survey  geoid. 

The  latitude  and  longitude  errors  from  the  true  position  were  first  determined. 

These  were  then  investigated,  by  partition,  using  the  TACTICS  command  ANALYZE.  This 
gives  a print  out  of:-  the  number  of  cases  examined,  the  minimum,  maximum,  mean, 
standard  deviation  and  standard  error. 


The  following  formulae  are  used  to  calculate  these  figures.  In  all  cases  n is  the 
number  of  readings,  x is  the  mean  value  of  the  readings  and  x^  is  the  j^*^  reading. 


Standard  Deviation  a = 


r 

I 


Standard  Error  = this  is  sometimes  called  the  Standard 

v'  n 

Deviation  of  the  means. 

The  results  can  be  found  in  Section  9. 

The  HISTOGRAM  command  was  then  used  to  produce  histograms  of  latitude  and  longitude 
for  each  partition.  These  can  be  seen  in  Appendix  C. 

4.1  Editing  by  the  Investigator 

On  examination  of  the  histograms,  it  was  apparent  that  there  were  several 
fixes  which  had  errors  much  larger  than  the  others.  There  were  four  such  records. 
Numbers  17,  18,  229  and  284. 

The  full  record  of  each  was  checked,  and  in  each  of  the  first  three  cases  it  was 
found  that  the  difference  in  offset  frequency  from  the  previous  record  was  greater 
than  three,  which  causes  errors.  In  the  last  case  it  was  noted  that  the  ship's 
speed  had  been  increased  from  zero,  the  correct  figure,  to  20  knots.  This 
error  was  attributed  to  this  large,  erroneous,  change  in  speed,  and  the  fix 
rejected. 

These  four  records  were  excluded  from  all  further  calculations. 

5.  FURTHER  ANALYSIS  AND  HISTOGRAMS 

From  the  results  obtained  from  this  and  previous  analysis,  it  was  decided  to 
investigate  the  records  further  in  an  attempt  to  determine  possible  trends  and 
relationships  in  the  error  patterns. 

To  facilitate  this  analysis  a program,  SAT  2,  was  written  and  is  shown  in 
Appendix  D.  This  read  data  from  SATDAT,  transformed  it  to  a useable  form  and 
deposited  in  a file  SATR. 

The  error  in  latitude  and  longitude  was  calculated,  in  nautical  miles.  This  used 
the  fact  that  one  minute  of  latitude  equals,  for  a first  order,  one  nautical  mile 
and  that  longitude  error  in  nautical  miles  equals  longitude  error  in  minutes 
multiplied  by  the  cosine  of  the  latitude. 

The  number  of  good  doppler  counts  before  and  after  the  centre  of  pass  were  also 
calculated. 

Finally  the  radial  error  was  calculated  from  the  equation:- 


/ . 2 2 
Radial  error  = / (latitude  error)  + (longitude  error) 

and  the  bearing  of  the  satellite  fix  position  from  true  North  was  calculated. 

The  TACTICS  package  was  now  brought  into  use  and  the  data  partitioned,  as  described 
in  Section  4. 
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The  data  was  printed  out  by  partition  and  is  contained  in  Appendix  E. 

An  analysis  and  histogram,  by  partition,  was  produced  for  latitude  error,  longitude 
error,  radial  error  and  bearing  from  North.  The  analysis  can  be  found  in  Section  9 
and  the  histograms  in  Appendix  F. 

5.1  Histogram  Distributors 

The  histograms  were  then  investigated,  and  in  some  cases  further  histograms 
with  smaller  ranges  were  plotted  so  that  their  shape  could  be  seen  more  clearly. 

These  histograms  can  be  found  in  Appendix  G. 

5.1.1  Latitude  Error  and  Longitude  Error  Histograms 

The  laticude  and  longitude  histograms  for  the  updated  fixes  (partition 
FIX. UPDATE)  show  clearly  defined  peaks  and  the  points  distributed  in  a 
normal  manner.  The  mean  was  not  zero  thus  showing  there  is  an  offset 
from  the  actual  latitude  and  longitude. 

The  latitude  and  longitude  histograms  for  the  fixes  not  used  for  a 
position  update  (partition  FIX. NOUPDATE) , and  Velocity  North  fixes 
(partition  VNFIX),  cover  a wider  range  due  to  the  greater  scatter  of 
the  points.  Because  of  this  the  majority  of  the  points  tend  to  be 
within  three  or  four  intervals.  Again  the  means  were  not  at  zero,  thus 
indicating  a position  offset. 

In  order  to  investigate  further,  more  histograms  were  plotted.  These 
were  centred  on  the  mean  and  had  a range  of  about  one  standard  deviation. 
These  histograms  can  be  found  in  Appendix  G. 

To  confirm  that  the  distributions  were  Gaussian  (normal),  a normalised 
Gaussian  distributor  was  calculated  from  the  measured  mean  and  standard 
distribution.  A Chi-squared  test  was  then  carried  out  between  the 
measured  cumulative  distribution  and  the  calculated  values.  The 
Chi-squared  figure  was  always  sufficiently  small  to  be  able  to  say  that 
the  measured  distribution  was  part  of  the  calculated  distribution.  In 
all  cases  points  beyond  three  standard  deviations  from  the  mean  were 
ignored.  This  is  statistically  valid  as  these  points  are  probably  due 
to  systematic  and  not  random  errors. 

5.1.2  Radial  Error  and  Bearing  from  True  Position  Hi stograms 

With  a perfect  radial  error  distribution,  affected  only  by  random 
noise,  a half  normal  distribution  would  be  expected,  with  a peak  at 
zero. 

The  histogram  for  the  updated  fixes  shows  a half  normal  distribution 
but  the  peak  is  offset  to  a distance  of  about  0.035  nautical  miles. 

This  offset  was  to  be  expected  from  the  mean  latitude  and  longitude 
cal culations . 

The  histograms  for  the  two  other  partitions,  the  NO-UPDATE  fixes  and 
the  Velocity  North  fixes,  shows  a large  peak  and  a large  range,  being 
5 nm  and  4 nm  respectively.  When  the  large  peak  was  plotted  in  greater 
detail  a histogram  of  the  same  shape  as  that  for  the  updated  fixes  was 
produced. 

The  single  frequency  satellite  navigation  system  is  less  precise  in  its 
longitude  fixes  than  its  latitude  fixes,  and  will  give  a wider 
distribution  of  fixes  in  the  East-West  direction  than  in  the  Nortn-South 
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direction.  Thus  a histogram  of  bearing  from  North  would  be  expected 
to  show  two  peaks,  one  at  90°  (corresponding  to  East)  and  one  at  270° 
(correspondi nq  to  West).  It  must  be  remembered  that  on  the  bearing 
histogram  360°  is  equivalent  to  0°.  The  bearing  histogram  for  updated 
fixes  did  not  conform  to  that  predicted.  Only  one  peak  was  apparent 
and  this  was  centred  at  a bearing  of  about  280°. 

The  N-UPDATE  fix  histogram  was  more  scattered,  probably  due  to  the 
small  number  of  cases  used,  but  a peak  could  be  seen  at  about  300°. 

The  Velocity  North  histogram,  on  the  other  hand,  showed  the  same 
distribution  as  the  updated  fixes,  with  only  one  peak  centred  at  about 
290°.  This  clearly  called  for  further  investigation,  and  is  dealt  with 
in  the  following  section. 

6.  SCATTER  PLOTS 

A scatter  plot  is  a plot  of  latitude  against  longitude;  from  this  one  can 
observe  the  dispersion  of  the  fixes.  Plots  for  all  three  partitions  were  produced, 
and  these  can  be  found  in  Appendix  H;  this  appendix  also  contains  the  program  output 
for  each  partition.  A plot  routine,  SAT  3,  was  written  to  produce  the  scatter  plots. 
This  program  requires  the  operator  to  type  in:  a title,  the  outer  limit  of  the  plot, 

the  file  holding  the  input  data  and  the  file  to  which  the  plot  data  is  to  be  written. 
During  operation  the  program  types  out  any  position  fixes  which  are  exactly 
coincident  and  any  fixes  which  are  outside  the  outer  limit. 

The  program  used  subroutines  from  Comshare's  CIL  Plotter  package,  and  produced  data 
suitable  for  driving  a CIL  plotter.  Since  the  Admiralty  Compass  Observatory  does  not 
have  a plotter  capable  of  generating  the  scatter  plots,  paper  tapes  were  generated 
by  the  Comshare  computer  and  taken  to  ASWE  where  a suitable  device  is  held.  The 
plotting  routine  is  listed  in  Appendix  I. 

The  plot  of  updated  fixes  demonstrated  that  most  of  them  lie  to  the  North  and  West 
of  the  true  position,  thus  accounting  for  the  bearing  of  280°.  It  could  then  be 
seen  that  the  dispersion  was  greater  in  the  East-West  direction  than  in  the  North- 
South  di'-ection.  Thus  the  explanation  of  the  bearing  histogram  is  that  the 
displacement  of  the  mean  from  the  true  position  is  sufficiently  great  to  conceal  the 
East-West  distribution. 

The  plot  for  the  N-UPDATE  fixes  shows  a more  scattered  distribution,  but  there  is 
still  a noticeable  bias  to  the  West  and  North,  and  the  East-West  distribution  is 
larger  than  the  North-South  distribution. 

The  Velocity  North  partition  scatter  plot  shows  very  clearly  the  North  and  West 
bias  with  a larger  East-West  distribution  than  North-South  distribution. 

Thus  the  North  and  West  bias  shown  on  the  scatter  plots  explains  why  the  bearing 
histograms  had  peaks  at  about  290°. 

7.  STATISTICAL  ANALYSIS 


The  normal  distribution  of  the  latitude  and  longitude  histograms,  would  tend 
to  indicate  that  errors  are  due  to  a random  (white  noise)  effect.  To  test  the 
hypothesis,  the  TACTICS  package  was  used  to  correlate  all  the  variables.  Three 
correlation  matrices  were  produced,  one  being  generated  for  each  partition.  These 
matrices  are  contained  in  Appendix  J. 
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The  correlation  coefficient  is  equal  to:- 

n n n 
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where  n is  the  number  of  readings 

X is  the  latitude  or  longitude  reading 

y is  the  parameter  of  the  value  against  which  the 
correlation  coefficient  is  required. 


For  the  UPDATE  fixes,  the  highest  correlation  factor  occurred  when  comparing 
Latitude  errors  against  offset  frequency.  At  0.3207,  this  figure  is  not 
sufficiently  high  to  establish  any  significant  correlation  between  Latitude  error 
and  offset  frequency. 


The  other  two  partitions  show  significant  correlation  between  Latitude  error  and 
offset  frequency,  for  the  Velocity  North  fixes  there  is  also  a correlation  between 
Longitude  error  and  offset  frequency. 


To  investigate  this  further,  graphs  of  latitude  error  and  longitude  error  against 
offset  frequency,  for  the  particular  partitions  in  question,  were  produced, 
(Appendix  K). 


In  all  cases  it  can  be  seen  that  there  is  a cluster  of  points  around  the  most 
common  offset  frequencies,  but  that  any  large  offset  frequency  has  a corresponding 
large  error. 


Thus  the  correlation  does  not  denote  a linear  relationship  between  latitude  error 
or  longitude  error  and  offset  frequency,  but  simply  confirms  the  fact  that  a large 
offset  frequency  will  produce  a large  positional  error.  Having  explained  tin's 
there  are  no  other  significant  correlations,  thus  the  postulated  hypothesis  must 
have  a high  degree  of  probability. 

To  verify  this  beyond  doubt,  the  linear  stepwise  regression  program  contained  in 
the  TACTICS  package  was  used  on  both  latitude  errors  and  longitude  errors  by 
partition.  They  were  regressed  onto  DAY,  Minutes  (MIN),  Iterations  (ITER), 
Elevation  (ELEV),  0-ffset  frequency  (FREQ),  number  of  doppler  counts  before  centre 
(NOBC)  and  number  of  doppler  counts  after  centre  (NOAC).  The  term  R-Square  is  a 
measure  of  fit:  the  closer  it  is  to  unity  the  better  the  fit.  This  is  defined  as 

the  proportion  of  the  total  variation  about  the  mean  of  the  dependant  variable, 
which  is  accounted  for  by  the  regression  equation. 


The  results.  Appendix  L,  show  that  only  the  Velocity  North  fixes  have  any  significant 
value  of  R-squared,  and  that  for  longitude  the  error  was  only  marginally  significant. 
For  both  latitude  error  and  longitude  error,  the  major  part  of  the  fit  was  due  to  the 
offset  frequency. 

It  is  therefore  possible  that  there  is  some  linear  relationship  between  Velocity 
North  position  errors  and  offset  frequency.  It  is  still  not  possible  to  say  that 
there  is  a definite  relationship  as  over  10%  of  the  distribution  about  the  mean 
remained  unaccounted  for.  This  is  not  good  enough  to  make  a definite  statement 
about  a relationship. 
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8.  INVESTIGATION  OF  EQUIPMENT  PRODUCED  STATISTICS 


This  equipment  produces  three  positional  fix  accuracy  indicators:- 

S-LA  - the  standard  deviation  of  the  latitude  in  minutes  of  arc 
based  on  the  rms  Doppler  residuals. 

S-LO  - the  standard  deviation  of  the  longitude  in  minutes  of  arc 
based  on  the  rms  Doppler  residuals. 

RMS  - the  calculated  rms  error  of  the  fix,  in  nautical  miles  using  the 
previous  two  quantities  and  the  quantity  S-VN  (the  calculated 
standard  deviation  of  the  Velocity  North)  where  applicable. 

In  doing  this,  the  equipment  manufacturers  have  made  the  assumption  that  the 
majority  of  the  latitude  and  longitude  errors  can  be  attributed  to  the  Doppler 
calculation  residuals.  In  order  to  test  if  this  was  a valid  assumption,  the 
values  of  the  measured  latitude  and  longitude  errors,  in  minutes  of  arc,  and  the 
absolute  values  were  calculated  and  then  correlated,  by  partition,  with  the 
respective  three  equipment  values  (Appendix  M). 

For  updated  fixes,  there  is  no  significant  correlation  between  the  measured  values 
and  those  calculated  by  the  equipment.  For  the  two  other  partitions  a fairly 
high  degree  of  correlation  is  shown  between  the  measured  values  and  equipment 
produced  figures. 

To  investigate  further,  graphs  were  produced  of  measured  absolute  latitude  and 
longitude  errors  against  the  equipment  calculated  values  (Appendix  N). 

On  investigation  each  graph  showed  a similar  pattern  consisting  of  most  of  the 
points  clustered  near  the  origin  but  with  a small  number  of  points  a great  distance 
away.  The  conclusion  that  may  be  drawn  from  this,  is  that  very  large  residuals, 
which  result  in  the  equipment  producing  a large  value  of  S-LAT  and  S-LON,  do  give 
bad  fixes  and  thus  large  latitude  and  longitude  fix  errors.  This  is  a fact  which 
is  already  known. 

It  was  now  necessary  to  investigate  if  there  was  any  linear  relationship  between 
measured  and  equipment  produced  values.  To  do  this  a further  TACTICS  facility, 
regression,  was  used.  This  allows  a multivariate,  stepwise,  linear  regression  to 
be  performed  on  any  variable,  the  resulting  outputs  are  in  Appendix  0.  R-squared 
is  a measure  of  how  well  the  measured  points  fit  the  equation,  the  closer  it  is  to 
unity  the  better  the  fit. 

For  the  update  fix  latitude  and  longitude  errors,  R-square  is  very  low.  This  agrees 
with  the  correlation  figures  and  indicates  no  linear  relationship  between  measured 
and  equipment  produced  figures. 

The  other  partitions  show  a significant  value  for  R-squared,  except  for  the 
longitude  errors  in  the  N-UPDATE  fixes  where  R-square  has  a medium  to  low 
significance  of  0.37. 

For  N-UPDATE  fixes,  the  measured  latitude  and  longitude  errors  are  equal  to  some 
constant  plus  about  three  times  the  corresponding  equipment  calculated  values. 

Thus  for  this  type  of  fix,  the  equipment  tends  to  underestimate  the  error. 

For  Velocity  North  fixes,  the  measured  latitude  and  longitude  errors  are  equal  to 
some  constant  plus  about  one  third  the  correspondi ng  equipment  calculated  values. 

In  this  case  the  equipment  appears  to  have  over  estimated.  In  truth  the  equipment 
has  assumed  that  the  Northerly  velocity  is  unknown,  and  has  put  it  as  a variable 
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in  its  calculations.  Thus  although  the  set  was  stationary,  it  has  calculated  a 
Northerly  velocity  and  therefore  produced  a larger  calculated  latitude  and 
longitude  error  (S-LA,  S-LO). 

It  is  therefore  apparent,  that  the  equipment  calculated  "error  indicators"  are 
markedly  different  from  the  measured  figures.  As  an  indication  of  fix  quality, 
they  could  therefore  be  quite  misleading  and  of  little  use  to  the  navigator. 

9.  SUMMARY  OF  RESULTS 

9.1  Latitude 


Actual 

Latitude  51°  29. 

45'  North. 

9.1.1 

Partition  Fix; 

UPDATE 

Minimum  51° 

29.39' 

(Error  = -0.06 

nm) 

Maximum  51° 

29.586' 

(Error  = 0.136 

nm) 

Mean  51° 

29.476' 

(Error  = 0.0250 

nm) 

Standard  Deviation  0.0211' 

(=  0.0211  nm) 

Standard  Error 

0.00178' 

(=  0.00178  nm) 

No  of  readings 

140 

9.1 .2 

Partition  Fix: 

N-UPDATE 

Minimum  51° 

29.31 ' 

(Error  = -0.14 

nm) 

Maximum  51° 

29.582’ 

(Error  = 1.132 

nm) 

Mean  51° 

29.519' 

(Error  = 0.0687 

nm) 

Standard  Deviation  0.204' 

(=  0.204  nm) 

Standard  Error 

0.0286' 

(=  0.0286  nm) 

No  of  readings 

51 

9.1.3 

Partition  VN  Fix 

Minimum  51° 

28.238' 

(Error  = -1.212 

nm) 

Maximum  51° 

30.502' 

(Error  = 1 .052 

nm) 

Mean  51° 

29.485' 

(Error  = 0.0355 

nm) 

Standard  Deviation  0.151' 

(=  0.151  nm) 

Standard  Error 

0.0137' 

(=  0.0137  nm) 

No  of  readings 

121 

) 
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Actual 

Longitude  0°  33, 

.583'  West. 

9.2.1 

Partition  Fix: 

UPDATE 

Minimum  0° 

33.474' 

(Error  = 

-0.0679  nm) 

Maximum  0° 

34.01  ' 

(Error  = 

0.266 

nm) 

Mean  0° 

33.701 ' 

(Error  = 

0.0735 

nm) 

Standard  Deviation  0.0874' 

(=  0.0544 

t nm) 

Standard  Error 

0.00738' 

(=  0.0046 

i nm) 

No  of  readings 

140 

9.2.2 

Partition  Fix; 

N-UPDATE 

Minimum  0° 

30.372' 

(Error  = 

-1 .999 

nm) 

Maximum  0° 

41.078' 

(Error  = 

4.642 

nm) 

Mean  0° 

33.761 ' 

(Error  = 

0.111 

nm) 

Standard  Deviation  1.239' 

(=  0.771 

m) 

Standard  Error 

0.173' 

(=  0.108 

nm) 

No  of  readings 

51 

9.2.3 

Parti tion  Fix: 

VN  Fix 

Minimum  0° 

29.286' 

(Error  = 

-2.676 

nm) 

Maximum  0° 

39.162' 

(Error  = 

3.474 

nm) 

Mean  0° 

33.697' 

(Error  = 

0.071 

nm) 

Standard  Deviation  0.752' 

(=  0.468 

nm) 

Standard  Error 

0.064' 

(=  0.0426  nm) 

No  of  readings 

121 

4 
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Radial 

Error  from  True 

Position  1 

9.3.1 

Partition  Fix: 

UPDATE  1 

Mi nimum 

0.0207  nm 

Maximum 

0.297  nm 

Mean 

0.0846  nm 

Standard  Deviation 

0.0479  nm 

Standard  Error 

0.00405  nm 

No  of  readings 

140 

9.3.2 

Partition  Fix: 

N-l 

UPDATE 

Mi  nimum 

0.0218  nm 

Maximum 

4 . 7 1 8 nm 

Mean 

0.346  nm 

Standard  Deviation 

0.729  nm 

Standard  Error 

0.102  nm 

No  of  readings 

51  y 

9.3.3 

Partition  Fix: 

VN 

FIX  I 

Mi nimum 
Maximum 
Mean 

Standard  Deviation 
Standard  Error 
No  of  readings 


0.0282  nm 
3.679  nm 
0.201  nm 
0.456  nm 
0.0414  nm 
121 


13 


9.4  Angular  Displacement  from  North 
9.4.1  Partition  Fix:  UPDATE 


Mi  nimum 

o 

o 

o 

Maximum 

348.72° 

Mean 

276.67° 

Standard  Deviation 

65.879° 

Standard  Error 

5.5678° 

No  of  readings 

140 

Partition  Fix: 

N-UPDATE 

Minimum 

10.081° 

Maximum 

350.92° 

Mean 

236.72° 

Standard  Deviation 

97.607° 

Standard  Error 

13.668° 

No  of  readings 

51 

Partition  Fix: 

VN 

FIX 

Mi nimum 

1.6971° 

Maximum 

352.02° 

Mean 

264.13° 

Standard  Deviation 

86.99° 

Standard  Error 

7 . 9082° 

No  of  readings 

121 
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10.  CONCLUSIONS 

The  system  performed  well  giving  good  accurate  fixes.  The  position  errors 
were  due  to  both  an  offset  error,  which  was  to  the  North  and  West  of  the  actual 
position,  and  a normally  distributed  random  variation. 

The  offset  error  could  be  due  to  the  actual  position  not  being  known  sufficiently 
accurately,  or  a difference  in  the  geoidal  reference  system  used  by  the  satellite 
navigator  computer  program  and  the  British  Ordnance  Survey  reference  system.  This 
is  being  investigated  with  the  Ordnance  Survey  Authorities. 

For  good  fixes  in  both  latitude  and  longitude,  the  offset  error  is  greater  than  the 
standard  deviation  of  the  noise.  A significant  increase  in  accuracy  could  therefore 
be  obtained  if  the  offset  error  could  be  removed. 

A statistical  analysis  of  the  results  demonstrated  that  there  was  no  linear 
correlation  between  the  accuracy  of  the  fix  and  any  of  the  variables  printed  out  by 
the  system.  It  was  observed  that  large  shifts  in  some  of  the  variables,  such  as 
offset  frequency,  did  produce  large  errors;  an  effect  which  was  postulated. 


The  investigation  of  the  equipments  estimate  of  standard  deviation  for  latitude  and 
longitude,  based  on  the  position  calculation  residuals,  showed  that  in  the  case  of 
good  fixes  there  was  no  relationship  between  these  figures  and  the  actual  measured 
figures.  There  was  some  degree  r.  correlation  for  bad  (N-UPDATE)  fixes  and  Velocity 
North  fixes,  and  this  varied  from  about  three  times  too  large  for  Velocity  North 
fixes  to  three  times  too  small  for  bad  fixes.  In  consequence  a navigator  could  put 
little  reliance  on  the  equipment  produced  fix  parameters  and  as  such  they  are  very 
misleading. 

The  Velocity  North  fixes  were  on  the  whole  worse  than  the  standard  fixes.  This 
would  indicate  that  adding  Velocity  North  as  an  unknown  to  the  program,  has  resulted 
in  an  increase  rather  than  a reduction,  in  the  calculated  position  error.  It  must 
be  remembered  that  all  results  were  obtained  with  the  receiver  stationary,  and 
therefore  no  vehicular  motion  effects  were  introduced.  The  purpose  of  the  VN  fix 
is  to  improve  the  VN  assessment  of  the  vehicle  motion,  and  in  a dynamic  situation 
where  there  may  well  be  SPEED  LOG  errors  some  improvement  may  well  exist. 

11.  RECOMMENDATIONS 


The  equipment  produ  .u  standard  deviation  and  rms  figures  are  largely 
erroneous,  and  thus  misleading.  These  figures  should  be  removed  from  the  printout. 
They  could  be  replaced  by  a fix  quality  figure  based  on  the  residuals  of  the 
computer  calculation.  Navigators  would  then  learn  from  experience  the  magnitude  of 
the  errors  associated  with  each  fix  quality  figure. 

Magnavcx  should  thoroughly  investigate  the  Velocity  North  program.  Since  the 
manufacturer  has  not  released  details  of  the  VN  program,  no  meaningful  suggestions 
can  be  made.  Nevertheless,  if  the  program  does  treat  Velocity  North  as  a complete 
unknown,  taking  no  account  of  the  input  speed  (either  automatic  or  manual),  this 
could  be  a problem  area.  Before  final  publication  of  this  report  Magnavox  has 
studied  its  contents  and  a Manufacturer's  Note  concerning  this  recommendation  has 
been  included.  (See  Page  16). 

The  equipment  supplier  should  be  asked  to  provide  the  geoidal  co-ordinate  system 
used  in  the  computations.  With  this  information  it  would  be  possible  to  calculate 
the  theoretical  offsets. 

Twelve  equipments  are  now  in  service  with  the  Royal  Navy,  and  at-sea  experience  is 
now  being  gained.  In  order  to  supplement  the  static  assessment,  the  in-service 
operators  should  keep  good  records  of  fixes  for  subsequent  analysis.  The  value  of 
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the  VN  solution  could  then  be  determined  in  a dynamic  environment.  It  is  possible 
that  this  may  be  simulated  by  deliberately  applying  faulty  velocity  inputs  into  a 
static  receiver,  and  this  should  be  further  studied.  Equipment  availability 
prevented  any  further  extension  of  the  assessment  described  in  this  report. 

12.  REFERENCES 

1.  P A C Kennedy  and  T A Jones.  "Performance  of  Satellite  Positioning  Corp 
Type  SCS-100  Satellite  Navigation  Equipment" 

ASWE  TR-74-19  (UNLIMITED). 


MANUFACTURER'S  NOTE 


The  following  comments  on  this  report  have  been  received  from  the  Magnavox  Co:- 

The  report  is  considered  to  be  good,  giving  a fair  assessment  of  the  operation  of 
Magnavox  systems  software.  Magnavox  would  like  to  point  out  that  the  comment  made 
in  the  report  as  to  the  Velocity  North  portion  of  the  programme  could  have  been  more 
meaningful  had  the  system  been  tested  by  entering  false  velocities  in  the  north  or 
south  direction.  This  would  cause  the  software  to  attempt  to  resolve  true  velocity 
north  and  give  a more  accurate  position  fix  as  a result. 

The  V North  programme  has  been  updated  to  take  care  of  problems  associated  with  the 
velocity  solution  and  steps  are  being  taken  to  retrofit  existing  systems  with  the 
corrected  software. 

The  report  comments  on  the  sigma  values  printed  out  after  the  fix  calculation  are 
being  considered  by  Magnavox  but  the  company  points  out  that  these  values  are  not 
provided  to  give  a firm  indication  of  fix  accuracy  but  only  as  a rough  guide. 
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Data  had  some  errors  but  the  computer  was  able  to  edit  them  out 


4i'h;i  FI  : HCQ-zv-rT'OSu 


hhv  TI^1E  LBT  ton  hdg  spd 

1 ?4  ll^eOO  051  i"?.4GS  N 0 00  5 3.74c  M 9 0.00  95.  00  . 0000 

I TEP  ElEV  GEQM  IhT  5-LH  S-LQ  5- Vf-I  RMS 

0 5 Z-;3.  Z-E  50190  .0199  . 0909  . 0000  . 0099 


Z-iJF'DHTE  -.01999790 
L“ i'lHE'o  i_— t'l.F’D 
17  0.  0 . ij  1 1 5 

iJPDrtTE 

01  . 0 9 5 9 9 4 04  . 0 0 0 0 4 0 9 'z 

09  - . 0 0 0 9 9 '9  75  - . i j 0 0 0 0 0 9 Z 

0 5 .0  0 9'  Z:  0 175  . 0 0 0 0 0 1 7 Z 


. 17999778 


0 0 0 01779  - . 0 0 0 0 0 4 3 9 
0 0 0 0 0 1 9 0 . 0 0 0 0 0 0 0 9 
0 0 0 0 LI  549  - . 0 0 0 0 0 0 1 9 


mill  0 0 0 09 


40  0 VMFI:;  hCD-Z:V-75030 


riH','  TIME  LHT  LGM  HMT  HDl-.  5PD 

1 Z'4  1 159.00  051  99.500  M 000  33.  393  W 90.00  93.  00  . OOL 

I TEC’  ELEV  GEDM  ShT  3-LR  3-LD  Z-Vh  RMS 

03  33.  S-E  30190  .0359  .0949  .9197  .0050 


Z-ijPDRTE  . 0159  0590 

i_— I'lHEiG  L-  — l-'IZF'Ei 


04  1 7 05  37 


. 04490009 


M-UPEiHTE 


Figure  1:  SECTION  OF  ORIGINAL  TELETYPE  LOG 


poor  fix. 
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APihNllIX  A LISTING  OF  PlAOGaAM/SAT''  AND  THE  PRIOTCUP  Pl-KiDUCSD  BY  I? 

• ^4l4lTAJ  / 3AT1-SY;..AT.NAV  THURSDAY  APR  24,  1975  9:00:59  AM 
'4  CARD 

C SATl  SORTS  THK  ROUGH  DATA  FZ<OK  FILE  SATUAT  IN'TO  FILE  SaTD 
REAL  LATM,  LONM 

INTEGER  DAY,  TLME,  SAT,  UPD,  DNS,  DEW,  TLAT,  TLON 
DI.NENSION  AUF(3),  TFIX(3) 

DaTA  AUP  /'UPDA  -UP  ' / , TLAT,  TLON  / 51 , 0 / 

D„TA  TFIX  / 'FIX  UNFI  ' / 

C REAjj  in  a record,  exit  'TO  1000  IF  NO  MORE  DaTA 
OPEN  (1,  INPUT,  / SATDAT-DA) 

OPEN  (2,  OUTPUT,  / SATD-DA) 

PRINT  ('SaTNAV  DATA  SORT  ROUTINE') 

ICNT=0 

10  READ  (1,  5,  END=100)  II,  TFIX(3) 

5 iXIRMAT  ( / 13,  IX,  A4) 

READ  (1,  15)  day,  tike,  LATD,  LATM,  DNS,  LOND,  LONM, 

& DEW,  aNT,  HDG,  SPD,  NA’EL 
15  FORMAT  ( / / / 31,  F,  A1,  I,  F,  A1 , 4F) 

KEaD  (1,20)  ITER,  ELEV,  GEOM,  SAT,  SLA,  oLO,  3VN,  i^\S,  FRE , 
20  FORMAT  ( / I,  F,  A3,  I,  5F) 

REaD  (1,25)  rtUP(3) 

25  FORf'.aT  ( / / / A4) 

C TEoT  THaT  LAT  aNU  LON^  ARE  CORRECT 

IF  (LATD.  S^.  TLAT.  AND.  LOND.  E«i.  TLON)  GuTO  27 
PRINT  ('LAT/LON  ERROR  ON  RECORD  = ' , I),  ICNT 
PRINT  ('LaT  = ',  I,  'LON=  ',  I),  LaTD,  LOND 

GOTO  35 


A1 


1 

1 

i 

C I’KoT  i’HnT  i’lX  13  OKAY 

27  Ir’(TrIX(.3).r:.,.rFIX(2).  OK.  rFIX(3)  .TFIXCl ) ) OOTO  28 

; PRINT  ('  KELP  TEIX,  HECOKU  = I),  ICNT 

i STOP 

I 

C TEST  THAT  UPD..TE  IS  OKAY 

28  IF(,iUP(5).E.^.AUP(2).LK.nUP(5).EH.-'.UP(l))  SoTO  29 

PRINT  (’iiliP  UPU,  RECOHL  = I),  ICNT 

STOP 

C SET  UP  OUTPUT  RECORD 

29  aRITE  (2,  50)  D«Y,  TIME,  LaTD,  LATM,  DNS,  L^ND,  LONM,  DEW, 

ANT,  HDG,  SPD,  NVEL,  TPIX(3),  AIJP(5), 

^ INTER,  ELEV,  OEOM,  SAT,  SLA,  SLO,  SVN,  RKS,  FRE.^ 

50  FOR^UT  (5(1,  IX),  F,  IX,  A1,  IX,  I,  IX,  F,  IX,  A1 , 4(1X,  F) , 

& 2(1X,  A4),  IX,  I,  IX,  F,  IX,  A5,  IX,  I,  5(1X,  F) ) 

55  ICNT  = ICNT  + 1 

SOTO  10 

100  PRINT  ('  NUMBER  OF  RECORDS  = I),  ICNT 

PRINT  ('  RECORDS  STORED  IN  / SATD-DAT') 

CLOSE 
STOP 
END 


A2 

i _ 


PHOGrUK  OUTPUT 


SATNAV  UAPA  oOKT  ROUTINE 
L,.T/LON  EPPOR  UU’  i-ffiCORD  = 19 
LAT  = 55  LON  = 6 
LAT/LON  ERROR  ON  iiECORD  = 20 
LAT  = 55  LON  = 8 
NUMBER  OF  RECOPLiS  = 3lS 
RECORDS  STORED  IN  /SmTD-D,>T 
(5  MnIN)  100S  + 5 


ERRONEOUS  RECORDS  IN  / SAT  DAT 

400  FIX  hCO-SY -75050 

000001 1111111111111111111111111111 100000 

A^ 

D«Y  TIME  LnT  LCK  ANT  HDG  SPD  NVEL 

065  004800  055  05.979  N 006  07.551  W 200.0.0000.0000.0000 
ITER  ELEV  SEOM  SAT  S-LA  S-LO  S-VN  RMS  FRE.^ 

04  29.  S-W  30159  28.52  46.14  .0000  16.93  - 3187.6 
S-UPDATE  216.91627  - 203.07232  300.57744 
C— V/HDG  C-w’SPD 

270.2  39758 

N-UPDATE 


400  VNFIX  ACO-SY-75030 

00001 1111111111111111111111111111 00000 

AA 

DAY  TIME  LAT  LON  ANT  HDG  NVEL 


063  004800  055  20.241  N 008  39.002  W 200.0.0000.0000  - 402.0 


02  29.  S-W  30130  176.6  350.0  397.9  18.49  -3044.6 
S-UPDaTt  231.2CA12  - 302.54082  380.77074 

C-iVHDG  G— rtSPD 

270.1  64821 


N-UFDATK 


APPEND  I X_  B:  Listing  of  Ddta  File  SAID 

62  1533C0  bl  29*544  n 0 33*629  w 2q0  0 0 0 FIX  N**uF  <= » 1 N-W  3014Q  *02?8  *07?2  0 *0018  - 1 0 * *»  1 
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/^..T-4Y;  ^,.TNnV  Fi<IUAY  JUN  2?,  1975  9:3^:  44  AM 

/{  CAA.J 

C /^,.T2  oOrfTo  THi;  .iOUGi!  DAT/.  FKOM  FILF  /ShTD.xT-Ijh 

C ItJTO  FILE  /S.xTH-Dn  GUGH  THAT  IT  IS  SUITABLE  FOB 

G ilEGRESSIO.x  WOHK, 

C 

ILixiiL  LATM,  LONM 

IL'TEGEH  DAY,  TIME,  ShT,  TLAT,  TLON,  HRS,  SATNUM(6) 

DI.'.ENoION  DOPPLER  (40),  CENTRE  (40) 

D.-xTA  0;E  /'I'/,  ARRO'w  /'/<'/,  TLAT,  TLON  /51,0/ 

DATa  TLaTH,  TLOxNM  /29. 450, 33.583/,  SERO  /'O'/ 

DxiTA  SaTNLTM  / 30120,  30130,  30140,  3OI8O,  30190,  30200/ 

d.ita  pi  /3. 141592654/ 

C OPEN  INPUT  .xND  OUTr’UT  FILES 

OI’EN  (1,  INPUT,  /.„xTDAT-D.x) 

OPEN  (2,  OUTPUT,  /S.xTxR-DA) 

C iRINT  P.ROaR,xM  NAME 

xRl,,T  ( / 'xROSxRxIM  /SAT2'  / 2X,  ' oORTS  D/xTa  FROM  /S.xTDaT-DA' , 
& 'INTO  .iEG.-iKoSION  D-xTa  FILE  /o,xTR-D-0 
C SET  NUMBER  OF  RECOiL^S  GOUNTEx!  TO  SEi^ 

ICKT=0 

C CLEAR  DOPi  LER  AND  CENTRE  ARRAYS 

10  DO  11  1=1,40 
DOPPLER  (I)  = EERO 

11  GENTxffi  (I)  = EERO 

c Rxi:a./  a record  of  twelve  li..e:; 

R.LxJ  TYPE  OF  FIX 
KE,x.-  (1,100,  END=1000)  Ti’''IX 

D1 


C Hiliij  DOPPLEH  DATA 

READ  (1,110)  DOPPLER 
C REhj  CENTRE  DATA 

READ  (1,110)  CENTRE 
C READ  DAY,  TIME  LAT  LONG 

READ  (1,120)  DAY,  Tlf-E,  LAT,  LATM,  LON,  LONM 
C READ  ITER  ELEY  SAT  & FREQ 

READ  (1,150)  ITER,  ELEY,  SAT,  XI,  X2,  X5,  X^t,  FREQ 
C READ  UPDATE 

READ  (I,l40)  UPDATE 

C TEST  THAT  LAT  AND  LONG  DEGREES  ARE  CORRECT 

IF  (LAT.  NE.  TLAT.  OR.  LON.  NE.  TLON)  GOTO  10 
C CALCUL.ATE  LAT  AND  LONG  ERRORS  IN  NAUTICAL  MILES 
ELAT  = LATM-rTLATM 

ELON  = (LONM-TLONM)  • COS ( (FLOAT (TLAT) +TLATM/6O. ) *(PI/l80.)) 
G CALCULATE  NUMBER  OF  MINUTES  AFTER  MIDNIGHT 
HRS  = TIME/10000 
MINS=TIME/1 00-HRS ♦ 1 00+HRS  * 60 
C CALCULATE  NUMBER  OF  DOPPLERS  BEFORE  AND  AFTER 
C CENTRE  OF  PASS 

C CLEAR  CENTRE  ARROW  REACHED  FLAG  NARR 
C CLEAR  NUMBER  OF  DOPPLERS  BEFORE  CENTRE  COUNTER  NOBC 

C CLEAR  NUMBER  OF  DOPPLERS  AFTER  CENTRE  COUNTER  NOAC 

NARR=0 
N0BC=0 
N0AC=0 

C START  CALCULATIONS 
DO  20  1=1,  40 


D2 


IF  (DOPPLEH  (I).  NE.  ONE)  GOTO  15 
IF  (NAM.  E^.  0)  N0BC=N0BC+1 

IF  (NARR.  Ey.  l)  N0AC=N0AC+1 

15  IF  (CENTRE  (I).  Ey.  ARROW)  NARR=1 
20  CONTINUE 

C CALCULATE  WHICH  SATELLITE 

DO  5 1=1,  6 

IF  (SAT.  E(^.  SATNUM  (I))  SAT  = I 
30  CONTINUE 

C CALCULATE  THE  FADIAL  ERROR  DISTANCE 

R=SyRT  (ELAT ‘ELAT  + ELON  *ELON) 

C CALCULATE  ANGLE  OF  POSITION  FROM  DUE  NORTH 

C REMEMBER  THAT  LONGITUDE  IS  WESTERLY 

ELON— ELON 
ANG=  0. 

IF  (ELAT. GT.O. 001 . OR  ELAT. LT.-O. 001 ) GOTO  40 
IF  (ELON. GT.O. 001)  ANG=90. 

IF  (ELON. LT. 0.001)  ANG=270, 

GOTO  70 

40  IF  (ELON. GT.O. 001.  OR. ELON. LT.-O. 001 ) GOTO  50 
IF  (ELAT.  LT.O.)  ANG=l80. 

GOTO  70 

C CALCULATE  ATAN  ANGLE 

50  TANG=ATAN(ELON/ELAT)  •(l80/Pl) 

IF(ELON.LT.O. ) GOTO  60 
ANG=TANG 

IF(ELAT.LT.O. ) ANG=l80.+  TANG 
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BEARING  FROM  TRUE  NORTH  HISTOGRAM,  FIX  : N- UPDATE  PARTITION 


I 

i 


0 

18 

f\> 

13 

> 

36 

o 

36 

> 

54 

o 

54 

72 

o* 

72 

> 

90 

cn 

90 

108 

— 

108 

> 

126 

O 

126 

> 

144 

o 

144 

> 

162 

o 

lo2 

> 

180 

o 

130 

> 

198 

o 

133 

216 

— 

216 

> 

234 

23't 

> 

252 

ro 

232 

> 

270 

oi 

270 

> 

288 

u> 

288 

> 

306 

<£> 

306 

324 

U* 

324 

> 

342 

2 36 


BEARING  FROM  TRUE  NORTH  HISTOGRAM,  PARTITION  V^FIX 


0 

18 

cn 

1C 

> 

36 

— 

3u 

54 

54 

> 

72 

04 

72 

90 

fj* 

90 

> 

CO 

o 

ro 

108 

> 

126 

o 

126 

144 

o 

144 

162 

o 

lo- 

;Co 

o 

180 

l.i8 

o 

198  > 216 


o 


APPENDIX  6 
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G SAT3*  THIS  li-iOGRAM  RE^OS  Li-iT  ERRORS  AND  DONS  ErRuRS 
C EROM  a user  ^PECIEIED  P'ILE  and  produces  a D.-iTA  eide 
c vjhich  can  be  used  by  the  cil  plotter 

G 

INTEGER  TITLE(3) 

DIME:JSI0N  EL.-Vr(350),  EL0N(350),  IEILE(5),  CEILE(5) 

C DEFINE  input  mND  OUTPUT  FILES 

PRINT  CFROGR-iM  SAT3'/,  3^,  'SCATTER  PLOT  ROUTINE'/) 
0PEN(2,  INPUT,  PROMPT  'INPUT  D.VTA  i’lLE  ',  lEILE) 
OPENd,  BIN;G(Y  output,  prompt  'PLOT  EI1.E  ',  UEILE) 
DEC0JE(i6, (5A4) , OEILE)  OEILE 
DEC0DE(i6, (5A4),  liTLE)  lEILE 

ChLL  plots  (1) 

C i<E.vD  IN  a TITLE 

DISPLnY  'PLOT  TITLE  ' 

AGCEi’T  110,  TITLE 

C iiEAD  INEIXEo  AND  COUNT  THE  NUMBER  READ 
ICNT=0 

10  ;-iEAD  (2,  100,  END=20)  ELAT (ICNT+1 ) , ELON(ICNT+l) 

ICNT=ICNT+1 

C LONGITUDES  jvRE  j^EST,  GET  THE  CORRECT  ORIENThTION 
ELON ( ICNT ) =-EL0N ( ICNT) 

SOTO  10 

20  DISPLAY  'OUTER  LIKIT(NK)  ' 

■ iCCEtT  100,  OL 

C T*^^»T  i*\)R  C0I.1' JIDEnCE  OF  f’CiNTS 


II 


JO  ho  1=1,  io;,T-i 
jO  30  J=1+1,  ICNT 

Ij?'(2Lm  ,l).;;i^.i:.L<v?(J}.C«.hLOrJ(l)..N£.£LON(J))  JO'i'O  30 

i-tiiNT( i’lii  a^Ms;  ',f.3x, 'lon  ',f,3x, 

& '!=',  1,3X,  ' iiL-vT(l),ELON(I),I,J 

30  cor.'TiNUi, 

40  OONTINUF 

C jRrtW  COMPAv,.a  LiriSo  rtiJiJ  Oli^CLFS 
O.-^i^  PLOTdO.,  10. ,-3) 

C-d.L  P0L(0.,0.,8.,4,/(,0) 

G .v;iOTATi;  TliE  GKID 

Gall  GYMBOL(0.,9.,0,5, 'NCHTH' ,0.,5) 

G/iLL  LUM3KH  ( 0 , , 8 . 25 , 0 . 5 , OL , 0 . , 2 ) 

C i’LOT  THE  POIKTG 

DO  60  I=1,1GNT 

i C T^jOT  IF  FOILT  10  .vITHIN  THE  LIMITS 

I IF(S..HT(ELaT(I)**2+El0N(I)**2).LE.0L)  GOTO  50 

PHIHTCoUT  OF  i<ANGE  LnT  ',F,3X,'lON  ' ,F,3X, 

& 'I=',D,  ELAT(I),  ELON(I),  I 

GOTO  60 

50  X=(ELUiHI)/0L)*8 

Y=(ELAT(I)/0L)'8 

Call  ghaf(x,y,o.i .2) 

j oO  GOL'TxiJUE 

C PiIIaT  the  .'ITLE 

1 C<vLL  GYKB0L(-3.,-10.,0.5,TITLE,0.,12) 

*'  G FIL'IGH  PLOTTING 

G,»LL  ILOT  (0.,0.,999) 
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Appendix  L:  Output  from  stepwise  regression  program  to  test  position  errors  against 

day,  minutes,  iterations,  elevation,  offset  frequency,  number  of  Doppler  counts 
before  centre  and  number  of  Doppler  counts  after  centre  of  pass. 


> 


REG>  STE  EUT  ON:  VI -V7  BYP 


LATITUDE  ERROR 


PARTITION  FIX. UPDATE 


EQUATION 

EUT  = - . 07073  + .00178*DAY  + .00490*FREQ 

R-SQUARED  = .13874 


PROCEDURE  SUMMARY 


VARIABLE 

R-SQUARED 

F-RATIO 

STEP 

ENTERED  REMOVED 

VALUE 

% CHANGE 

CRITERION 

ACTUAL 

1 

FREQ 

.1028 

3.91 

15.82 

-> 

DAY 

. 1387 

34 . 90 

3.91 

5.71 

PROCEDURE  ENDED  WHEN 

MAXIMUM  F TO  ENTER  = .63  FOR  MINS 

FAILED  TO  MEET  F ACCEPT  = 3.91 


LATITUDE  ERROR 


PARTITION  FIX. NOUPDATE 


EQUATION 

EUT  = - 1.6329  + .02179*DAY  - .05514*FREQ 

R-SQUARED  = .56864 


PROCEDURE  SUMMARY 

VARIABLE  R-SQUARED  F-RATIO 

ACTUAL 

50.95 
6.55 

PROCEDURE  ENDED  WHEN 

MAXIMUM  F TO  ENTER  = 1.63  FOR  NOAC 

FAILED  TO  MEET  F ACCEPT  = 4.05 


STEP 

ENTERED 

REMOVED  VALUE 

CHANGE 

CRITERION 

1 

FREQ 

.5098 

4.04 

2 

DAY 

.5686 

11.55 

4.04 

-LI- 


LATITUDE  ERROR 


I’ARTITION  V'NFIX 


EQUATION 

ELAT  = - .56871  + .00727*DAY  - .02907*FREQ 

R-SQUAREl)  = .88802 


PROCEDURE  SUMMARY 


VARIABLE 

R- 

SQUARED 

F -RATIO 

STEP 

ENTERED  REMOVED 

VALUE 

“b  CILXNGE 

CRITERION 

ACTUAL 

1 

FREQ 

.8760 

3.92 

840.59 

DAT 

.8880 

1.37 

3.92 

12.68 

PROCEDURE  ENDED  MIEN 

MAXIMUM  F TO  ENTER  = 1.84  FOR  ELEY 
FAILED  TO  MEET  F ACCEPT  = 3.92 


REF>  STE  ELON  ON;  V1-V7  BYP 


LONGITUDE  ERROR 


PARTITION  FIX. UPDATE 


PROCEDURE  SUNMARY 

NO  INDEPENDENT  VARIABLE  QUALIFIED  FOR  ENTRY 

PROCEDURE  ENDED  WHEN 

MAXIMUM  F TO  ENTER  = 2.30  FOR  FREQ 

FAILED  TO  MEET  F ACCEPT  = 3.91 


LONGITUDE  ERROR 


PARTITION  FIX. NOUPDATE 


PROCEDURE  SUMMARY 

NO  INDEPENDENT  VARIABLE  QUALIFIED  FOR  ENTRY 

proci;dure  ended  when 

MAXIMUM  F TO  FJMTER  = .62  FOR  ELEY 

I'AILF.D  TO  MEET  F ACCEPT  = 4.04 


-L2- 


LONtUTUDK  tRROR 


PARTITION  VNFIX 


EQUATION 

ELON  = .18602  + .16387*ITER  + .07753*FREQ  - .01791*N0BC 

R-SQUARED  = .68966 


PROCEDURE  SUMMARY 


VARIABLE  R-SQUARED  F-RATIO 


STEP 

ENTERED 

REMOVED 

VALUE 

% CHANGE 

CRITERION 

ACTUAL 

1 

FREQ 

.6612 

3.92 

232.27 

7 

NOBC 

.6731 

1.79 

3.92 

4.27 

3 

ITER 

.6897 

2.46 

3.92 

6.25 

PROCEDURE  ENDED  WHEN 

MAXIMUM  F TO  ENTER  = .55  FOR  DAY 

FAILED  TO  MEET  F ACCEPT  = 3.92 


.VPPINDIX  M:  CORRELATION  MATRIX.  TOH  EACH  PARTmON,  OP  MEASURED  POSITION 

ERRORS  EITH  EQUIPMENT  CALCULATED  ERRORS 

PARTITION  FIX.UPDATE 
140  OBSERVATIONS 


LATE 

LONE 

ALAT 

ALON 

SLA 

SLO  PMS  fl 

LATE 

1 .0000 

lone 

.1430 

1 .0000 

J 

ALAT 

.831? 

.1342 

1 .0000 

s 

ALON 

.1223 

.9642 

.1030 

1.0000 

1 

SLA 

- .0102 

.01433 

.2021 

- .0202 

1 .0000 

1 

SLO 

- .0347 

.0284 

-.0482 

.0260 

- .1021 

1.0000  1 

RMS 

.0351 

- .0277 

.0481 

.0253 

.1006 

-1.0000  1.0000  1 

f 

PARTITION  PIX.NOUPDATE 
51  OBSERVATIONS 


NOTE: 


LATl'RJDE  ERROR  IN  DEGREES 

LONGITUDE  ERROR  IN  DEGREES 

ABSOLUTE  VALUE  OP  LATITUDE  ERROR 

ABSOLUTE  VALUE  OP  LONGITUDE  ERROR 

EQUIPMENT  CALCULATED  LATITUDE  ERROR 

EQUIPMENT  CALCULATED  LCMGITUDE  ERROR 

EQUIPMENT  CALaiLATED  FMS  POSITION  ERROR  IN  NM 


M - 1 


LATE 

LCNE 

ALAT 

ALON 

SLA 

SLO 

RMS  ii 

LATE 

1 .0000 

1 

LONE 

.5353 

1 .0000 

ALAT 

.9513 

• 5175 

1 .0000 

1 

ALON 

.4947 

.5654 

.5306 

1 .0000 

Q 

SLA 

.8866 

.4543 

.9513 

.4572 

1 .0000 

n 

SLO 

.1838 

- .1842 

.2565 

.6052 

.2566 

1 .0000 

j 

RMS 

• 5861 

.8147 

.6031 

.8943 

.5354 

.2797 

1.0000  9 

PARTITION  UNPIX 

1 

121  OBSERVATIONS 

LATE 

LONE 

ALAT 

ALON 

SLA 

SLO 

RMS  1 

LATE 

1 .0000 

1 

LONE 

- .8449 

1 .0000 

ALAT 

-.1375 

.1536 

1 .0000 

1 

ALON 

- .2072 

.0617 

.8658 

1 .0000 

1 

SLA 

.2879 

- .2590 

.8588 

.7587 

1 .0000 

1 

SLO 

.2179 

- .4298 

.6274 

.7935 

.7589 

1 .0000 

1 

RMS 

.2765 

- .2460 

.0904 

.0872 

.1152 

.1955 

1 .0000  1 
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' APPENDIX  0:  OUTPUT  PROM  SfEPWISE  REGRESSION  PROGRAM  TO  TEST 

MEASURED  ATiP  EQUIPMENT  CALCULATED  POSITION  ERRORS 
ABSOLUTE  LATITUDE  ERRORS 
PARTITION  FIX.  UPDATE 
EQUATION 

^ ALAT  = .02254  + .45623*SLA 

; R-SQUARED  = .04086 

i 

PROCEDURE  SUMMARY 


STEP 

ENTERED 

VARIABLE 

REMOVED 

R-SQUARED 
VALUE  io  CHANGE 

CRITERION 

F-RATIO 

ACTUAL 

1 

SLA 

.0409 

3.91 

5.88 

PROCEDURE  ENDED 

WHEN 

MAXIMUM  P TO  ENTER  = .11  FOR  RMS 
FAILED  TO  MEET  P ACCEPT  = 3.91 

PARTITION  FIX.  NOUPDATE 

EQUATION 

ALAT  = -.01750  + 2.5i68*SLA  + .80966*RMS 
R-SQUARED  = -91730 

PROCEDURE  SUMMARY 


VARIABLE 

R-SQUARED 

P-RATIO 

STEP 

ENTERED 

REMOVED 

VALUE 

$ CHANGE 

CRITERICH 

ACTUAL 

1 

SLA 

.9050 

4.04 

466.6' 

2 

RMS 

.9173 

1.36 

4.04 

7-1' 

PROCEDURE  ENDED  WHEN 
MAXIMUM  P TO  ENTER  = .02  FOR  SLO 
FAILED  TO  MEET  F ACCEPT  =4.05 


0 - 1 


PARTITION  VJIPX 


EQUATION 

ALAT  = .0^245  + .39674*SLA 

R-SQUARED  = .73757 

PROCEDURE  SUMMARY 


VARIABLE 

R-SQUARED 

P-RATIO 

STEP 

ENTERED  REMOVED 

VALUE  fo  CHANGE 

CRITERICN 

aciua: 

1 

SLA 

•7376 

3.92 

334.41 

PROCEDURE  ENDED  WHM 
MAXIMUM  P TO  ENTER  = .63  FOR  SLO 

FAILED  TO  MEET  P ACCEPT  = 3.92 


ABSOLUTE  LONGITUDE  ERROR 


PARTITION  PIX.UPDATE 
PROCEDURE  SUMMARY 

NO  INDEPENDENT  VARIABLE  QUALIFIED  FOR  ENTRY 

PROCEDURE  ENDED  WHEN 
MAXIMUM  P TO  ENTER  = .09  for  SLO 
FAILED  TO  MEET  P ACCEPT  = 3.91 

PARTITION  PIX.NOUPDATE 

EQUATION 

ALON  = -.01102  - 1.5203*SLA  + 1.933l*SLO  + 30.450*RMS 
R-SQUARED  = .94212 

PROCEDURE  SUMMARY 


VARIABLE  R-SQUARED  P-RATIO 


STEP 

ENTERED 

REMOVED 

VALUE  $ 

CHANGE 

CRITERION 

ACTUAL 

1 

RMS 

• 7997 

4.04 

195.62 

2 

SLO 

• 9365 

17^11 

4.04 

103.42 

3 

SLA 

• 9421 

.60 

4.05 

4.56 

PROCEDURE  ENDED  WHEN 
ALL  ORGINAL  VARIABLES  E.'TERED 


0-2 


PARTITION  VNPIX 


EQUATION 

ALCN  = .07143  + 2.3665*ALAT  + .21135*LCXJE  + .89591 *SLA 

R-SQUARED  = .53853 


PROCEHJRE  SUMMARY 


STEP 

VARIABLE 

ENTERED  REMOVED 

R-SQUARED 
VALUE  $ CHANGE 

P-RATIO 

CRITERION 

] 

1 

1 

ACTUAL  ’ 

1 

1 

ALAT 

.5024 

3.87 

313.03 

2 

LONE 

.5140  2.30 

3.87 

7.34 

3 

SLA 

.5385  4.78 

3.87 

16.38 

PROCEDURE  ENDED  WHEN 

MAXIMUM  P TO  ENTER  = 1.97 

FOR  LATE 

FAILED  TO  MEET  P ACCEPT  =3.87 

[I 
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0 - 3 
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